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Abstract 



This paper complements the study of single top production at the LHC aiming to 
estimate the sensitivity of different observables to the magnitude of the effective couplings. 
In a previous paper the dominant VF-gluon fusion mechanism was considered, while here 
. . ■ we extend the analysis to the subdominant (10% with our set of experimental cuts) s- 

; I ' channel process. In order to distinguish left from right effective couplings it is required 

^ to consider polarized cross-sections and/or include nii, effects. The spin of the top is 

accessible only indirectly by measuring the angular distribution of its decay products. 
We show that the presence of effective right-handed couplings implies necessarily that 
the top is not in a pure spin state. We discuss to what extent quantum interference terms 
can be neglected in the measurement and therefore simply multiply production and decay 
probabilities clasically. The coarsening involved in the measurement process makes this 
possible. We determine for each process the optimal spin basis where theoretical errors 
are minimized and, finally, discuss the sensitivity in the s-channel to the effective right- 
handed coupling. The results presented here are all analytical and include mt, corrections. 
They are derived within the narrow width approximation for the top. 
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1 Introduction 



At present not a lot is known about the Wtb effective coupling. This is perhaps best evidenced 
by the fact that the current experimental results for the (left-handed) Ktb matrix element 
give in 

\K 

' 0.99 ±0.29. (1) 



\Ktd\^ + \Kt,\^ + \Ktb\^ 

In the Standard Model this matrix element is expected to be close to 1. It should be empha- 
sized that these are the 'measured' or 'effective' values of the CKM matrix elements, and that 
they do not necessarily correspond, even in the Standard Model, to the entries of a unitary 
matrix on account of the presence of radiative corrections. These deviations with respect to 
unitary are expected to be small — at the few per cent level at most — unless new physics 
is present and makes an unexpectedly large contribution. At the Tevatron the left-handed 
couplings are expected to be eventually measured with a 5% accuracy [Q]. 

As far as experimental bounds for the right handed effective couplings is concerned, the 
more stringent ones come at present from the measurements on the b ^ decay at CLEO ^ . 
Due to a mt/rrih enhancement of the chirality flipping contribution, a particular combination 
of mixing angles and effective right-handed couplings can be bound very precisely. The 
authors of Q reach the conclusion that |Re((7ij,)| < 0.4 x 10~^. However, considering gji as 
a matrix in generation space, this bound only constraints the tb element. Other effective 
couplings involving the top remain virtually unrestricted from the data. The previous bound 
on the right-handed coupling is a very stringent one. It should be obvious that the LHC will 
not be able to compete with such a bound. Yet, the measurement will be a direct one, thus 
ruling out some contrived models where substantial cancellations might hypothetically avoid 
the 6 ^ S7 constraint. For the value of the effective couplings in some specific models see 
e.g. ||. 

At LHC energies the mechanism underlying single top production, therefore allowing a 
direct test of the Wbt effective couplings ql and gji, consists of several different processes 
(see e.g. [^). The dominant process is the so-called VF— gluon fusion channel, or t-channel 
process. The electroweak subprocesses corresponding to this channel are depicted in Fig. |^, 
where light u-type quarks or d-type antiquarks are extracted from the protons. Besides this 
dominant channel (250 pb at LHC [^) single tops are also produced through the process 
where the boson interacts with a 6-quark extracted from the sea of the proton (50 pb)[0] 
and in the quark-quark fusion or s-channel process (10 pb) which is depicted in Fig.|2[ The 
numbers quoted here correspond to total cross-sections. The separation between the sub- 
dominant processes and the dominant VK-gluon fusion is purely kinematical0, |^. By placing 
a cut on the pT of the detected b quark, the former process can be eliminated altogether. This 




Figure 1: Feynman diagrams contributing single top production subprocess. In this case we 
have a d as spectator quark 




Figure 2: Feynman diagram contributing to single top production in the subdominant s- 
channel process. The top decay is also shown in this figure 



also eliminates a sizeable fraction of the tops produced via the VF-gluon fusion mechanism 
(about two thirds for the cuts we use) . The cut on px has the additional bonus of making the 
QCD corrections manageable. One is therefore left with those single tops coming from the 
VF-gluon fusion mechanism (t-channel) and the subdominant s-channel process. The later 
one is actually the main object of our interest in this article, although we will also have many 
comments to make on the i-channel process. 

In a proton-proton collision a bottom-anti-top pair is also produced through analogous 
subprocesses. The analysis of such anti-top production processes is similar to the top ones 
and the corresponding cross sections can be easily derived doing the appropriate changes. 

In a previous paper |^ we have analyzed the sensitivity of different LHC observables to 
the magnitude of the charged current effective couplings considering only the dominant W- 
gluon fusion channel. In that work we did not consider the subsequent decay of the top 
in any detail. We did, however, a complete analytical calculation of the subprocess cross 
sections, for general left and right effective couplings and including bottom mass corrections. 
A > 30 GeV cut in the transverse momentum of the produced h quark was implemented in 

and, accordingly, only the so-called 2^3 process was retained, excluding top production 
off a 6-quark from the proton Fermi sea. Given the (presumed) smallness of the right handed 
couplings, the bottom mass plays a role which is more important than anticipated, as the 
mixed crossed gLdR term, which actually is the most sensitive one to qr, is accompanied by a 
h quark mass. The reader is encouraged to see where a very detailed analysis is presented. 

Typically the top quark decays weakly well before strong interactions become relevant, 
so we could in principle 'measure' its polarization state with virtually no contamination of 



strong interactions (see e.g. y, |10|] for discussions this point) and try to establish interesting 



observables based on this measurement. In fact it is not difficult to convince oneself that in 
order to disentangle left from right effective couplings, it is almost compulsory to be able to 
'measure' the polarization of the top. This will become apparent from the formulae presented 
in section ^. For this reason we have derived in this work and in[^ analytical expressions for 
the cross sections for the production of polarized tops or anti-tops. To this end one introduces 
the spin projector 

1 + 75 , 



with 

= —^^^=^^^= {pi ■ h,p\n) , = 1, = -1, (2) 

V (Pi)^ - {Pi ■ nf 

as the polarization projector for a particle or anti-particle of momentum pi with spin in the h 
direction. The calculation of the subprocess cross sections have been performed in this work 
and in P] for an arbitrary polarization vector n. 



Obviously, however, the top decays very shortly after production, so the only practical 
way one can measure the spin of the top is through its influence on the angular distribution 
of the leptons produced in the decay. It is tacitly assumed in most of the works published 
on this subject that the decaying top is in a pure spin state for all practical purposes; i.e. its 
polarization vector is pointing in a particular direction in space in a given reference frame. 

In the tree-level Standard Model this is not quite true, but it is almost true. The tree 
level Standard Model corresponds in our notation to taking gi = ^ and = 0. Imposing 
the a cut on px we have mentioned, only two subprocesses contribute; VF-gluon fusion and 
the s-channel process. The later provides 100% polarized tops in a certain direction (to be 
discussed latter). The situation in the t-channel process is a bit more complicated. The 
results from our previous analysis presented in |^ show that single top production is highly, 
but not fully, polarized in this case too (84 % in the optimal basis, with the present set of 
cuts). This is a high degree of polarization, but still well below the 90+ claimed by Mahlon 
and Parke in |jT^. We understand this being due to the presence of a 30 GeV cut in pT- In 
fact, if we remove this cut completely we get 91 % polarization, in rough agreement with ||l^ 
(note that we do not include the 2 — > 2 or b-sea process). Inasmuch as they can be compared 
our results for the tree-level Standard Model are in good agreement with those presented in 
1^ in what concerns the total cross-section. These considerations are quite independent of 
the choice of the strong subtraction scale, which is by far the largest source of uncertainty]^. 
Let us assume now for the sake of discussion that the polarization is indeed 100% . The 
top subsequently decays (say emitting a positively charged lepton). One can compute the 
angular probability distribution of the lepton with respect to the polarization direction in the 
Standard Model, multiply the two probabilities and compare the experimental result with 
the theoretical prediction. 

In fact things are a lot more subtle. First of all, we have seen that even in the Standard 

Model polarization is never 100% . Furthermore, it turns out that when g^ ^ 0, i.e. beyond 

the Standard Model, the top can never be 100% polarized (see the discussion in section Q 

and in ^]), not even in principle. In other words, the top is necessarily in a quantum mixed 

state and is described by a density matrix. The entries of this density matrix depend on the 

^ Since we perform a leading order calculation in QCD, the scale dependence is large. We have made two 
different choices: (a) fi — p™* is used as scale in as and the gluon PDF, while the virtuality of the W boson 
is used as scale for the PDF of the light quarks in the proton. This gives an excellent agreement with the 
calculations in [Q]. (b) jj,^ = s, s being the center-of-mass energy squared of the qg subprocess. The total cross 
section above the cut is then roughly speaking two thirds of the previous one, but no substantial change in 
the distributions takes place. This is the typical error for LO calculations in the present kinematical regime. 
The total cross section has been known to NLO for some time , while NLO results for the differential cross 
section have become available just recently jl2l 



momenta of the incoming and outgoing particles; that is to say, there is an entanglement 
between spin and momenta. 

Of course this complication amounts to a small effect because gR is surely quite small, 
even in most models beyond the Standard Model, so in first approximation the experimental 
consequences should be small. However, if our purpose is precisely to measure gji or at 
least to set a bound on it, it is clear that the effect needs to be taken into account. As 
already emphasized, to be able to tell left from right effective couplings one absolutely needs 
to consider the spin of the top. 

The next step is to select the direction where one is to 'measure' the spin of the top. 
By tracing the appropriate spin operator with the density matrix one would determine the 
expected probabilities of finding a top that (after the measure) would point in the given 
direction of our choice. There is a privileged spin basis, namely the one where the density 
matrix is diagonal, where the calculation is greatly simplified since one needs not compute 
the off-diagonal terms. This diagonalization process has to be done event by event and it 
selects a particular vector n (event dependent). In section || we provide explicit formulae for 
this privileged direction. Elementary Quantum Mechanical considerations show that this is 
also the direction where the differential cross section is maximal (or minimal depending on 
the sign of the spin). Using this 3- vector as spin basis, for instance, one can multiply the 
probability of producing a top polarized in the positive +n direction with the corresponding 
decay angular probability distribution plus the probability of producing a top polarized in the 
negative —n distribution times the the corresponding decay angular probability distribution. 
The dependence on the effective left and right couplings and gji is obviously contained in 
the density matrix and also in the decay distributions. 

Obviously, since the entries of the density matrix depend on the spin basis, the final 
physically observable result of the previous analysis will certainly depend on the spin basis 
too. How is this possible? In fact this is as it should be; we are multiplying probabilities 
and in fact we are neglecting the quantum interference terms because we are assuming that 
the polarization of the top is measured in the intermediate state before the top quark decays. 
Then there should be no surprise in that the interaction between the top and the apparatus 
measuring its spin modifies the final physical results. 

However, a proper measure of the top spin before it decays is impossible; the only way 
we learn about top polarization is precisely from the final decay products. So, the previous 
procedure it is conceptually incorrect^. The final result has to be strictly independent of the 
intermediate spin basis one uses. Does this mean that the usual procedure — which is the one 
we just described — is totally flawed? In principle yes, however one expects that the coarsening 
■^Even if one is considering, as we do here, only on-shell tops 



involved in the measuring process washes some or all of the interference effects. Then perhaps 
the previous procedure where one assumes that the spin of the top is well defined and one 
proceeds as if it could be measured before it decays it could be approximately correct. But 
what are then the errors involved? Do they jeopardize the determination of some of the 
effective couplings, in particular the distinction between gi and qr? These are some of the 
issues we would like to address in the present work. 



2 The differential cross section for polarized top production 

We shall discuss here the i-channel production for the sake of definiteness. This is the most 
involved process. We refer the reader to |^ for detailed expressions of the different amplitudes. 
We denote the matrix elements of the hard subprocess of Fig. || by M^. There will also be 
a M", corresponding to having instead a -u as spectator quark. We will also eventually 
define the matrix elements corresponding to the processes producing anti-tops as M", and 
M^. With these definitions the differential cross section for polarized tops da can be written 
schematically as 



da = P [ fu 



Mi 



+ fAMi\ 



where fu and denote the parton distribution functions corresponding to extracting a u- 
type quark and a d-type quark respectively and /? is a proportionality factor incorporating 
the kinematics and also the gluon distribution function. Using our analytical results for the 
matrix elements given in the appendix of [S] we obtain for the differential cross-section 



da = j3fu 



I |2 / , \ , I |2 /I, I L \ , 9*RgL + gR9*L, , ^ , ■aL9R-a*R3L, 

\9l\ {a + an) + \gR\ {b + bn) + — -{c + Cn)+i— ^—^ — dn 

I |2/ N , I |2/, , N , 9*r9l+9r9*l, ^ .9*l9r - 9*r9l , 
\gR\ [a - an) + \gL\ [b - bn) + — -[c-Cn)-t— — dn 



9l 9r 



A 



9L 

9R 



(3) 



where 



fu (a + an) + fd 



\fu {c -\- Cn -\- idn) -\- \ f(i {c — Cn — idn 



V i/h (c + Cn- idn) + \ fd {c-Cn+ idn) fu {b + bn) + fd{a- an) J ' 

(4) 

and where a, b, c, a„, bn, Cn and d„ are independent of the effective couplings gR and g^ 
and the subscripts n indicate linear dependence on the top spin four- vector n. All these 
quantities depend only on masses and momenta. The c, c„ and dn terms are proportional 
to the bottom mass and are therefore absent if one neglects m;, (this at first sight does not 



look unreasonable, given the energies involved). Inspection of the above differential cross- 
section reveals that in the nib limit, the only way to tell left from right effective couplings is 
precisely by considering and measuring polarized cross sections (the terms in a„, 6„) unless 
one is willing to rely strongly on the parton distribution functions^. For these reasons, both 
polarization and nib terms are quite important. 

We observe that A is an Hermitian matrix and therefore it is diagonalizable with real 
eigenvalues. Moreover, from the positivity of da we immediately arrive at the constraints 

det^ > 0, (5) 
TrA > 0, (6) 

that is 

(/„ (a + an) + fd(.b- hn)) {fu (b + 6„) + f^{a- an)) 

> \ {fu + + {cl + dl) {U - h)^ + 2cCn {fl - fl)) , (7) 

and 

ifu + fd) {a + h) + {fu - fs) {an + bn) > 0. (8) 

Note that it is not possible to saturate both constraints for the same configuration because 
this would imply a vanishing A which in turn would imply relations such as 

a + b _ /j- fu _ an - bn 
an + bn fd + fu a-b ' 

which evidently do not hold. Moreover, since constraints (0) and must be satisfied for 
any set of positive parton distribution functions we immediately obtain the bounds 



ab + anbn - ^ (c^ + + d^) > 



1 

anb + abn - -CCn 



b^ + a'-{bl + al) > i(c2-(c2+d2)). 

In order to have a 100% polarized top we need a spin four-vector n that saturates the 
constraint (|5|) (that is Eq.(|^) for each kinematical situation, that is we need A (n) to have a 
zero eigenvalue which is equivalent to have a unitary matrix C satisfying 

C^AC = diag (A, 0) , 

for some positive eigenvalue A. In general such n need not exist and, should it exist, is in any 
case independent of the effective couplings qr and ql- Moreover, provided this n exists there 



^The statement is exact if one uses the so-called effective W approximation, which is not terribly accurate 
for the present case and certainly not recommended [|l3|, but widely used in LIfC physics 



is only one solution (up to a global complex normalization factor q) for the pair {gn^gL) to 
the equation dcr = 0, This solution is just 



gi = aCu, 



gn = aC22- 



(9) 



Note that if one of the effective couplings vanishes we can take the other constant and 
arbitrary. However if both effective couplings are non-vanishing we would have a quotient 
gn/ gi that would depend in general on the kinematics. This is not possible so we can conclude 
that for a non- vanishing g^ ( g^ is evidently non- vanishing) it is not possible to have a pure 
spin state (or, else, only for fine tuned gn a 100% polarization is possible). 

Let us now give a very simple example to make the previous discussion more under- 
standable: in the un-physical situation where — > it can be shown that there exists two 
solutions to the saturated constraint (H), namely 



once we have found this result we plug it in the expression ^ and we find the solutions 
(0,5l) with gi arbitrary for the + sign and (5i?,0) with gn arbitrary for the — sign. That is, 
physically we have zero probability of producing a right handed top when we have only a left 
handed coupling and viceversa when we have only a right handed coupling. Note that in this 
case it is clear that having both effective couplings non-vanishing would imply the absence of 
100 % polarization in any spin basis. This can be understood in general remembering that 
the top particle forms in general an entangled state with the other particles of the process. 
Since we are tracing over the unknown spin degrees of freedom and over the flavors of the 
spectator quark we do not end up with a top in a pure polarized state. 

3 Cross-sections for top production and decay in the s-channel 

Let us now turn to the s-channel process. This is, as already mentioned, subdominant but 
non-negligible since it roughly amounts to a 10% of all single tops produced after our set of 
cuts are imposed. It is also a lot cleaner from a theoretical point of view, as QCD corrections 
are small. As a by-product we shall derive the differential decay width, which is applicable 
to both the t and s-channel processes. 

Using the momenta conventions of Fig. ^ and averaging over colors and spins of the initial 
fermions and summing over colors and spins of the final fermions (remember that we have 




(10) 



included a spin projector for the top) the squared amphtude for top production is given by 

o4Ar / 1 \2 



\Mr,. 



1 



9L 

9l9r + 9R9l 



2 1 Pi + mtn \ - X , I ,2 / Pi- , . , . 
9R\ ( qi 7, ) yl2 ■ P2) + \9l\ 92 7^ [qi ■ P2) 



[mt {qi ■ q2) + {q2 ■ Pi) {qi ■ n) - {q2 ■ n) {qi ■ pi)] 



, • 9L9r - 9R9l u a P a 

+ imi, i^- -ef,apan'^Piqiq2 



+ \9r\ 



2 / Pi + mtn \ ~ ^ , I ,2 

\9r\ [q2 — ] [qi ■ P2) + \9l\ 



pi - mtn . , 
qi 7. 1 (^2 • P2) 



j^^^9LgR + gR9L_ . ^^^j _^ . . _ . . ^^^j 



. 9l9r-9r9l ^ ^p^o. p a 



(11) 



where ^l, and are left and right couplings to light quarks and ql and gR are the effective 
couplings to the top- bottom system. In the numerical results we have taken gi = 1, gR = 
0; i.e. we stick to the tree-level Standard Model values in the light sector, but is quite 
straightforward to include more general couplings. Notice that, exactly as for the t-channel, 
the crossed gL9R terms vanish in the differential cross-section in the — > limit. Also for 
exactly the same reasons as in the t-channel analysis, modulo parton distribution functions 
effects, the differential unpolarized production cross section would be proportional to I^lP + 



\9r\ 



The differential cross section for producing polarized tops is then 



daf, 



f [xi,X2,iqi + q2f ,-^qcd) dxidx2 



1 



'4 jg^gi - 

d'pi dH2 ,2,, ,4e4, ^ , 



where / \ xi,X2, {qi + Q'2) :^QCDj dxidx2 accounts for the quarks parton distribution func- 
tions. 

The total decay rate of the top, on the other hand, with arbitrary left and right effective 
couplings is given by 



'VI/ 



|2 , I |2 

\9l\ +\9r\ 



9 9 9 (mr — m. 



-12mtmf, 



9L9*R + 9r9l 



mi + 



64TrmtPi 



The squared amplitude corresponding to the decay rate in the channel depicted in Fig. (||) 
summing over the top polarizations (with a spin projector inserted), averaging over its color 



and summing over colors and polarizations of decay products is given by 

|Mf I' = - A (gi ^ q2 ^ ki, p2 ^ -P2) , (12) 

where |M„p (gi k2, q2 ki, P2 — > — P2) is just expression (pl|) with the indicated changes 



in momenta. In the above expression g^, and are the left and right couplings corresponding 
to the lepton- neutrino vertex. We have assumed 5l = 1, 5_r = 0, but again this hypothesis 
can be relaxed. The decay rate differential distribution for this channel is given by 

\M^^f d'h d^k2 d'p2 ,4,4,, ^, ^ , 

dL n = K Q ^ 5 7: Q 7: (^TT Ki + Ko + P2 " Pi • 

2p1 (27r)=^ 2A;0 (27r)3 2A;0 (27r)3 2p0 ^ ' ^ 2 f2 i^ij 

Finally, using the narrow-width approximation, we have that the differential cross section da 
corresponding to Fig. ^ is given by 

da = ^danx'^. (13) 

±n 

4 The role of spin in the narrow-width approximation 

Within the narrow-width approximation we just discussed we decompose the process depicted 
in Fig. Q in two consecutive processes: the top production and its consecutive decay. In that 
set up we denote the single top production amplitude as A^^^^jj-p) and the top decay amplitude 
as i?p -|_^(p). In the polar representation we write 

A — \A p^f±(P^ 

B \ = \B ^-i ^l e*"^*^*"^ 

where p indicate external momenta and h {p) a given spin basis for the top. The differential 
cross section for the whole process is schematically given by 

da = I \A ,f^i \B ,f^i ) + A f^i )B f^^^^ (14) 



/I |2| |2 I \ |2| |2 

|^p,+n(p)| |^p,+n(p)| dp + |j4p^_^(p)| \Bp_^p-j\ dp 



Hence 



~l~2 y |^p,+n(p)| |^p,+n(p)| l^p,— ri(p)| |^p,— fi(p)| 

X cos (93+ (p) - If- (p) + L0+ (p) - CJ_ (p)) dp 

/I |2| |2 I i |2| |2 

|^p,+n{p)| |-Sp,+n(p)| dp + |ylp^_fj(p)| |i3p__^(p)| dp. (15) 

Since the axis with respect to which the spin basis is defined is completely arbitrary da is 
independent on this choice of basis. However within the narrow width approximation one 



never computes da following formula (0). The commonly used procedure |T^, 14] consists in 
computing the probability of producing a polarized top and then multiplying this probability 
by the probability of a given decay channel (see Eq. ([l3|)). This procedure is equivalent to 
the neglection of the interference term in formula ( [T5| ) as indicated there. First of all, as 
discussed in the introduction, if one neglects the interference term, the result depends on 
the spin basis; i.e. on the direction one chooses to measure the third component of the top 
spin. This is of course acceptable if one really performs a physical measure of the spin (in 
the n{p) direction in this case) since the interaction with the apparatus modifies the state. 
A dependence on the spin frame is however unacceptable if the spin is not measured before 
the top decays. 

Let us see whether this approximation can justified nevertheless. Clearly, the integration 
over momenta enhances the positive-definite terms in front of the interference oscillating one. 
If in addition we make a choice for h {p) that diagonalizes the top spin density matrix and 
thus maximizes |Ap _|_^(p-)| and minimizes |74p _^(-p)|, then we expect the interference term to 
be negligible when compared to / |^p,4.ft(p)| | dp even for small amount of phase 

space integration. In the s-channel we will see in the next section that the limit of ^ 
there exists a spin basis n (p) where |Ap is strictly zero. This basis is given by 



2 



n= — 7 rq2 - Pi 

mt V {l2 ■ Pi ) 



From this it follows that for small gji if we use that basis the interference integrand is 

I 2 I I 2 

already negligible with respect to the dominant term / | i?p _|_fj(p) | dp. For qr ^ 

one can still find a basis that maximizes |Ap ^_^(p)| (and minimizes |^p,-n(p)|) ^^'^ therefore 
diagonalizes the top density matrix p. In the next section we will show how to obtain such 
a basis that will be the one used in our numerical integration. In these simulations we have 
checked numerically that this basis is the one that maximizes da and therefore, on the same 
grounds, the one that minimizes the interference term. The same considerations can be 
applied to the t-channel process. 

Given that the observables are strictly independent of the choice of spin basis only if the 
interference term is included, we can easily assess the importance of the latter by checking 
to what extent a residual spin basis dependence is present. We have checked numerically 
this point by changing the definition of the spin basis n (p) and noting that our results are 
actually only weakly dependent on the choice of n (p) even for a small amount of coarsening. 
A 4% maximum variation in px distributions was found between the optimal diagonal basis 
and another basis orthogonal to the beam axis (that is, almost orthogonal to all momenta). 
Moreover we have checked that if spin is ignored altogether (by considering unpolarized 
top production) roughly the same amount of variation with respect to the diagonal basis is 



observed. Thus we conclude that even though the dependence on the choice of spin basis 
is not dramatic, its consideration is a must for a precise description using the narrow-width 
approximation taking into account the presumed smallness of the effective coupHng to be 
measured and how subtle the experimental distinction of left and right couplings turns out 
to be. 

5 The diagonal basis 

As stated in the previous section in order to calculate the top decay we have to find the 
basis where the polarized single top production cross section is maximal. The can do this 
maximizing in the 4-dimensional space generated by the components of n constrained by 

n ■ pi = 0, V? = —1, (16) 

where pi is the top four-moment, that is 

11 9 9 9 9 

Q n Pi + n PI + n PI 
Pi 



( 0\2 / o\2||^||2 / 1 1 , 2 2, 2 2\2 

KPi) = \V\) \\n\\ - [n p^+n pi+n p^) 



where ||n|| = y (n^)'^ + (n^)^ -|- (n^)"^, that is = ||n|| with n the normalized spin three- 
vector. From above equations we obtain 



n = 



{Pi) ~ ("-Vi + n'^Pi + h'^Pi)' 
n = n 



from which Eq. (^) follows immediately. Let us now find the polarization vector that maxi- 
mizes and minimizes the differential cross section of single top production. 

5.1 The £-channel 

We will begin with the i-channel the was analyzed in the previous chapter. Using Eq. (^ 
we define 

ttn = n ■ a, bn = n ■ b, 

Cn = n ■ c, dn = n ■ d, (17) 
and using Lagrange multipliers Ai and A2 for constraints (|l^) we maximize 

£7 + Ai (n^ + 1) + ■ Pi, 



obtaining the equations 



n 



2Ai 

A 

2Ai 



fu 

fd 



I ,2 ,1 ,2, , a*R9L+gRgl , .glgn- g*R9L , 

\9l\ a+\gR\ b + -^ i^c + t-^ — d 



\gRf a+\g,fb+ 9r9l+9r9l^ + ^9l9R- 9r9l ^ 



2A 



(18) 

(19) 
(20) 



and thus using Eqs. (|18| ) and 
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and therefore 
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with the normalization factor Ai given by Eq. (19). Note that in the idealized case fu 
f^ = f we obtain 

(a -b) -pi 



n = a 



mf 



-pi -{a-b)} , 



mt 



where a is the normalization constant that does not depend on / or the effective couplings. 
In the SM (gR = 0) we obtain 

n = a(^fu {-:^P^ -a^ - /j {^^-^Pi - b 
where a is a normalizing factor. 

5.2 The s-channel 

The s-channel differential cross section has the form 



da = P {fufd + fcfs) biP (as + an) + \9r\^ (bs + br, 
+ 3n9L^9R9l ^ ^ /l9R-9\9l^-^ 



where again /3 is a proportionality incorporating the kinematics, and where fu,c and fg^ 
denote the parton distribution functions corresponding to extracting a n, c-type quarks and a 



d, s-type quarks respectively. Using again the decomposition ( |17D and proceeding analogously 
to the i-channel calculation we obtain 

n = aUgil — —Pi-a] + \gR\ — ^Pi - o 



m 



+ ^^^^ t '^^^ ( : ^^'^ ( ^Pi -d]\, (21) 



where a is the normalizing factor that in this case (unlike in the t-channel result) does not 
depend on the parton distribution functions. From Eq. (^) we obtain 

af" = -mtqt^ {qi ■ P2) , 

= +mtqi {q2 ■ P2) , 

= +mb{q'^ {q2 ■ pi) - q2 {qi ■ Pi)) , 

= -mte^^^,pUlq2, 

hence replacing in Eq. ( pT] ) we arrive at 

= a^lgif {{qi-p2){q2-Pi)Pi - iqi-P2)mfq^2) 

+ \9r\^ [{<i2 ■P2) {qi ■pi)Pi - {q2 ■ P2) m^tii) 
9*r9l + gngl / u / n u /■ \\ 

+mbmt— {q'{ (92 ' Pi) - ^2 (^i ' Pi)) 

+ ^'-^^^^^^^mbm,e^ ^^^p^H } > (22) 

which is the basis we use in our numerical simulations. If we neglect qr we obtain 

_ I {qi ■ P2) {q2 ■ Pi) Pi - {qi ■ P2) mlqi^ 
n — it , 

Y {qi ■ P2f {q2 ■ Pif mf - {qi ■ p2f mfq^ 
where we have included the normalization factor and since Q'2 = the above reduces to 

, f 

mtn = ± -q2 - Pi 

\{q2-Pi) 



which is the result we have quoted in the previous section coinciding with ||10[] 



6 Numerical analysis of s-channel single top production 

Let us start this section by discussing the experimental cuts we have implemented. Due to 
geometrical detector constraints ||T5[| we cut off very low angles for the outgoing particles. The 
charged particles in the final state have to come out with an angle in between 10 and 170 
degrees to be detected. These angular cuts correspond to a cut in pseudorapidity \r]\ < 2.44. 



In order to be able to separate the jets corresponding to the outgoing particles we implement 
isolation cuts of 20 degrees between each other. These are the appropriate cuts for general 
purpose experiments such as ATLAS or CMS. 

The set of cuts used in this work are compatible with the ones used in the i-channel. Since 
in the previous paper |^] top decay was not considered, the equivalence is only approximate 
and a more detailed phenomenological analysis will be required in due course (it is actually 
quite straightforward with the help of the results presented here to redo the i-channel study, 
but this goes beyond the scope of this paper). The present analysis should however suffice in 
any case to identify the most promising observables and get a rough estimate of the precision 
that it can be reached. 

We use a lower cut of 20 GeV in the b jet^. This completely eliminates top production 
from a 6-quark from the proton sea and greatly reduces higher order QCD contributions. In 
the i-channel reduces the cross section to about one third of its total value, since typically the 
b quark comes out in the same direction as the incoming gluon and a large fraction of them 
do not pass the cut. Similarly, pT > 20 GeV cuts are set for the top and spectator quark jets. 
These cuts guarantee the validity of perturbation theory and will serve to separate from the 
overwhelming background of low pj- physics. These values come as a compromise to preserve 
a good signal, while suppressing unwanted contributions. They are very similar to the ones 



used in [0] and pO|. 

In order to calculate the cross section a of the process pp — > tb we have used the CTEQ4 



set of structure functions |17] to determine the probability of extracting a parton with a given 
fraction of momenta from the proton. To calculate the total event production corresponding 
to different observables we have used the integrating Monte Carlo program VEGAS We 
present results after one year (defined as 10^ seg.) run at full luminosity in one detector (100 
fb"^ at LHC). 

Since in order to be able to perform the effective Wtb coupling one definitely needs to tag 
the two 6-type quarks, this value for the luminosity is surely too high. The b physics program 
at ATLAS |jl^ , for instance, it is planned to be done at one tenth of the total luminosity to 
avoid pile-up effects. This is even more so in a dedicated detector such as LHCbp|[p^] where 
the appropriate figure is expected to be 2 fb~^. ATLAS plans to do most of the 6-physics runs 
before full luminosity is reached, for instance. We have nonetheless used the high luminosity 

*In the previous paper Q the value used was 30GeV. We have decided to use this lower value here to have 
a larger total cross-section without compromising the theoretical accuracy 

^This type of analysis is anyway not well suited for such a detector. The rapidity for LHCb is in the range 
1.6 < rj < 4.9 and the angular separation cut between jets imposed here is totally unfeasible. Furthermore, 
jet reconstruction is not possible. Clearly the implementation of this type of physics to this detector requires 
a lot more ingenuity. 



figure since at tliis stage tlie experimental strategy is not totally settled yet. 

The way we proceed is the following. We analyze the kinematics of each event including a 
b that passes the experimental cuts and reconstruct the vector n using the analytic formulae 
presented in the previous sections. As the reader will remember, this provides us with a 
spin basis that minimizes the quantum interference terms. We then proceed to multiply the 
probabilities classically — ^just as if we pretend that the top spin has been measured in the 
direction determined by h and we determine the decay probability distribution. We retain 
only those final states that pass the remaining cuts. 

In the same way and choosing arbitrary spin directions we are able to see how much the 
physical results depend on the interference term. We have found a 4% difference between the 
worst case (a spin direction perpendicular to almost all 3-momenta involved) and the optimal 
case (found analytically here). We have every reason to believe that, after the integration over 
momenta and the resulting coarsening, the interference term is in this basis all but negligible. 
The rest of the results presented in this section are all worked out in the optimal spin basis. 

Let us first review the results obtained in the framework of the tree-level standard model. 
This corresponds to taking gi = ^ and = in all our formulae. The results are summarized 
in Figs. ^, ^ |5|. The first of these figures show how the final products of the process are 
predominantly emitted in the axis direction (albeit not so much as in the case of the t-channel 
production) and in the same direction. The plot shows the direction respect to the beam of 
bottom and anti-bottom. Recall that a 10 degree cut is implemented, as well as a separation 
cut of 20 degrees among all jets. Fig. ^ shows the pT distribution for the 6, showing the 20 
GeV cut on the pT of the h enforced. 

Fig. P shows the invariant mass of the lepton and bottom system in the tree level Standard 
Model. Since we are working in the narrow width approximation, the distribution falls to zero 
just below the physical mass of the top and this refiects the part of the total momentum of the 
top carried away by the undetected neutrino. Figs. ^ and actually show the pT distribution 
for the bottom and lepton, respectively, that are produced in the top decay. As previously 
discussed, 20 GeV cuts on the respective pT are imposed. Even though some information is 
lost by the fact that the neutrino is not seen and therefore there is some amount of missing 
momentum, this does not seem to affect the sensitivity to the effective couplings too much. 
One could as well consider channels in which the ^ produced in the top decay, decays 
hadronically. In hadronic decays of the top a full reconstruction of the top mass would be 
feasible. 

Let us now move beyond the Standard Model. Since changing the value of ql (while 
keeping gj^ = 0) amounts to a simple rescaling, we shall concentrate on the more interesting 
case of varying gji. As a rough order-of- magnitude estimate for the effective gn coupling we 




Figure 3: Distribution of the cosines of the polar angles of the bottom and anti-bottom with 
respect to the beam line. The plot corresponds to single top production at the LHC with top 
decay included. The calculation was performed at the tree level in Standard Model. For the 
parton distribution functions we use /x^ = s = (gi + ^2)^- 
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Figure 4: Anti-bottom transversal momentum distribution corresponding to single top pro- 
duction at the LHC. The calculation has been performed at tree level in the SM {gi, = 1, 

gR = 0). 




Figure 5: Distribution of the invariant mass of the lepton (electron or muon) plus bottom 
system arising in top decay from single top production at the LHC. The calculation was 
performed at the tree level in Standard Model with fj? = s = {qi + ^2)^- 
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Figure 6: Bottom transversal momentum distribution corresponding to single top production 
at the LHC. The calculation has been performed at tree level in the SM (^l = 1, 5r = 0). 
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Figure 7: Lepton (electron or muon) transversal momentum distribution corresponding to 
single top production at the LHC. The calculation has been performed at tree level in the 
SM {gL = 1,9R = 0). 



take \gji\ = 5 x 10^^. This is still worse than the limit implied by 6 ^ 57, but is the sort of 
sensitivity that LHC will be able to set. The effects are linear in gj^, so it is easy to scale up 
or down the results. We have consider the possibility of qr having a phase and, accordingly, 
the experimental sensitivity to that phase. 

We have found that the anti-lepton plus bottom invariant mass distribution we just dis- 
cussed in the previous paragraph is actually sensitive to gR. Figs. ^ and ^ reflect this 
sensitivity with the second figure showing the statistical significance per bin. 




Figure 8: Event production difference between non-vanishing gji coupling caculations and the 
tree level SM ones {gn = 0). Differences are plotted versus the invariant mass of the lepton 
(electron or muon) plus bottom system arising in top decay from single top production at 
the LHC. We have taken gn = +5 x 10~^, +i5 x 10~^, —5 x 10~^ and —i5 x 10~^ in plots 
(a), (b), (c) and (d) respectively. With the present set of cuts, the total number of events in 
the Standard Model is 181,000. The total excess is 1,200, roughly 1% . 

We shall now show the dependence of the three pT distributions (b, b and lepton) to the 
modulus and phase of the effective coupling g^. In all cases the value g'^, = 1 is taken. The 
sensitivity to departures from the tree level SM is shown in Figs, (p^, ( pT| ) and ([l^). 



We also include the statistical significance per bin for the signal vs cos (Ofi) in Fig. (13) 
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Figure 9: Plots corresponding to differences (a), (b) (c) and (d) of Fig. divided by the 
square root of the event production per bin at LHC. The square of the quotient denominator 
can be obtained from Fig. (^) multiplying da/drriinv by the LHC 1-year full luminosity (100 
fb^^) and by the width of each bin (4 GeV. in Fig. (||)). Taking the modulus of the above 
plots we obtain the statistical significance of the corresponding signals per bin. Note that 
statistical significance has a strong and non-linear dependence both on the invariant mass 
and the right coupling gji^ However purely imaginary couplings are almost insensible to their 
sign. 




Figure 10: Statistical significance per bin with respect to anti-bottom transversal momentum. 
Like in Fig. @ we have taken gn = +5 x 10~^, +z5 x 10~^, —5 x 10~^ and — i5 x 10~^ in 
plots (a), (b), (c) and (d) respectively. Note that here statistical significance has a strong 
dependence on the anti-bottom transversal momentum but is almost linear on Re (gn) and 
almost insensible to the sign of Im{gR). 
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Figure 11: Like in Fig. (P) we have taken gn = +5 x 10^^, +i5 x 10"^, —5 x 10^^ and 
— i5 X 10^^ in plots (a), (b), (c) and (d) respectively. Note that here statistical significance 
has a strong dependence on the bottom transversal momentum and clearly favors positive 
values of Re (gji) and again is insensible to the sign of Im{g^. 




Figure 12: Statistical significance of the corresponding signals per bin with respect to lepton 
(electron or muon) transversal momentum. Like in Fig. (^) we have taken gn = +5 x 10~^, 
+i5 X 10^^, —5 X 10^^ and -ib x 10~^ in plots (a), (b), (c) and (d) respectively. Note that 
again statistical significance has a strong dependence on the lepton transversal momentum 
and clearly favors positive values of Re{gii) . The sign of Im{gii) cannot be distinguished. 



and vs cos {9tf,) in Fig. (14). cos {9ti) and cos {9tb) are the cosines of the angle between the best 
reconstruction of top momentum and the momenta of anti-Iepton and bottom, respectively. 
In these figures we can clearly see that low angles corresponds to bigger sensitivities. This is 
in qualitative accordance with Eq. (|12|) which, after inspection, tells us that anti-leptons are 
predominantly produced in the direction of the top spin and therefore most of those produced 
predominantly in the top direction come from a top mainly polarized in a positive helicity 
state. Thus the quantity of those anti-leptons is more sensitive to variations in gji. Even 
though this argument applies in the top rest frame, the fact that most of the kinematics lies 
in the beam direction makes it valid at least for this kinematics. With the cuts considered 
here, the Standard Model prediction at tree level for the total number of events at LHC with 
one year full luminosity is 180700. Using the values gL = 1, gR = +5 x 10^^ leads to an 
excess of 1220 events which corresponds to a 2.9 standard deviations signal. The gL = 1, 
gjl = —5 X 10~^ model has a deficit of 480 events which corresponds to a 1.1 standard 
deviations signal. Finally the 5l = 1, f?/? = x 10~^ model has an excess of 367 events 
which corresponds to a 0.86 standard deviations. We see that there is a large dependence on 
the phase of gR. 

It is perhaps interesting to remark that after considering the decay process, the sensitivity 
to gji is actually quite comparable to the one obtained in the t-channel, where it was assumed 
that the polarized top was observable. From this point of view, not much information gets 
diluted through the process of top decay. 

The implementation of careful selected cuts can slightly improve these statistical signifi- 
cances but since here we are interested in an order of magnitude estimate we will not enter 
into such analysis here. Moreover since backgrounds are bound to worsen the sensitivity the 
above results must be taken as order of magnitude estimates only. A more detailed analysis 
goes beyond the scope of this article. 



7 Conclusions 

In this paper we have performed a full analysis of the sensitivity of single top production in 
the s-channel to the presence of effective couplings in the effective electroweak theory. The 
analysis has been done in the context of the LHC experiments. We have implemented a set of 
cuts which is appropriate to general-purpose experiments such as ATLAS or CMS. The study 
complements the one presented in that was devoted to the dominant t-channel process. 

We have seen that the determination of the right effective coupling in such an experimental 
context is quite challenging. One has to include both polarization effects and corrections. 
Analytical formulae are presented. 




Figure 13: Statistical significance of the corresponding number of events per bin with respect 
to cos {9ti) = pi ■ {pi + Pb) I \pi\ \pi where pi and ph are respectively the tree momenta 
of the lepton (positron or anti-muon) and bottom. The combination pi + j% is the best 
experimental reconstruction of the top momentum provided the neutrino information is lost. 
Like in Fig. (^) we have taken ^/j = +5 x 10"^, +i5 x 10"^, — 5 x 10"^ and — i5 x 10"^ in 
plots (a), (b), (c) and (d) respectively. Note that again statistical significance has a strong 
dependence on cos {Qti)- 




Figure 14: Statistical significance per bin with respect to cos (Otb) = Pb'iPi + Pb) I \P[ \ \Pi + Ph\ 
where pi and ph are respectively the tree momenta of the lepton (positron or anti-muon) and 
bottom. The combination pi+Ph is the best experimental reconstruction of the top momentum 
provided the neutrino information is lost. Like in Fig. (^) we have taken = +5 x 10^^, 
+ih X 10~^, —5 X 10"^ and —ib x 10~^ in plots (a), (b), (c) and (d) respectively. Note that 
again statistical significance has a strong dependence on cos {Otb)- 



Unlike in the discussion concerning the single top production through the dominant t- 
channel, top decay has been considered. The only approximation involved is to consider the 
top 3jS Si real particle (narrow width approximation). 

We have paid careful attention to the issue of the top polarization. We have argued, 
first of all, why it is not unjustified to neglect the interference term and to proceed as if 
the top spin was determined at an intermediate stage. We have provided a spin basis where 
the interference term is minimized. A similar analysis applies to the i-channel process. We 
present here and explicit basis for this case too. We get a sensitivity to qr in the same 
ballpark as the one obtained in the i-channel (where decay was not considered). Finally we 
have obtained that observables most sensible to are those where anti-lepton and bottom 
momenta are cut to be almost collinear. 
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